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Abstract 

Periodic one dimensional hopping model is useful to study the motion of 
microscopic particles, which lie in thermal noise environment. The mean veloc- 
ity V/v and diffusion constant Djv of this model have been obtained by Bernard 
Derrida [J. Stat. Phys. 31 (1983) 433]. In this research, we will give the lim- 
its Vd and V/v and Dat as the number N of mechanochemical sates in 
one period tends to infinity by formal calculation. It is well known that the 
stochastic motion of microscopic particles also can be described by overdamped 
Langevin dynamics and Fokker-Planck equation. Up to now, the correspond- 
ing formulations of mean velocity and effective diffusion coefficient, Vl and 
in the framework of Langevin dynamics and Vp,Dp in the framework of 
Fokker-Planck equation, have also been known. In this research, we will find 
that the formulations Vd and Vl,Vp are theoretically equivalent, and numeri- 
cal comparison indicates that Dp),Dl, and Dp are almost the same. Through 
the discussion in this research, we also can know more about the relationship 
between the one dimensional hopping model and Fokker-Planck equation. 
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1 Introduction 

Many physical fT, '2] and biochemical phenomena, especially for the motion of 
motor protein [21 HI E], can be described by the periodic one dimensional hopping 
model. In this model, the particle jumps along a periodical linear track from one 
binding site to next one through the sequence of mechanochemical states [TJ [8]. 
The particle in state j can jump forward to state j + 1 with the rate Uj, or jump 
backward to state j — 1 with the rate Wj. After moving N sites forward the particle 
comes to the same mechanochemical state but shifted by a step size distance L (for 
motor protein kinesin L = 8.2nm). The dynamics of particle motion is described by 
the standard rate equations of occupation probabilities Pj{t) 



dp, it) 
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At steady state. 
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It's solution is 
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This model has been extensively studied [9l[T0] and its mean velocity and diffusion 
constant D]\[ has been obtained explicitly [TT] . 
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Besides one dimensional hopping model, the stochastic motion of particle also can 
be modelled by Langevin dynamics [121 [13] ^ind Fokker-Planck equation [T4l [8l [T5| [T6] . 
Intuitively, the Langevin dynamics and Fokker-Planck equation can be regarded as 
infinite mechanochemical states case of periodic one dimensional hopping model. In 
the framework of Langevin dynamics, the particle position is governed by the following 
Langevin equation 

where ^ is viscous friction coefficient, is Boltzmann's constant. = ^{x)—FextX, 
Fext is external load, $ is a (tilted) periodic potential with period L. T is absolute 
temperature and fsit) is Gaussian white noise. Based on the relation between the 
effective diffusion coefficient and the first two moments of the mean first passage 
time (MFPT), the mean velocity and effective diffusion coefficient can be expressed 
in quadratures ([17] [T8] ) 
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(9) 



= 0(L) — 0(0) and D is the free diffusion constant which satisfies the Einstein 
relation C,D = ksT. 
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At the same time, the stochastic motion of particle can be modelled by the Fokker- 
Planck equation [H US US] 



dp 
dt 
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dx 



C, dx dx 



0<x <L 



(10) 



where p{x, t) is the probability density for finding particles at position x and time t 
At steady state, 
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It's solution is (see [Tl[I7l[l9]) 

Jexp 
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is the probability fiux. So the mean velocity of particles is 
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In [20] we have known that the formulation of the effective diffusion coefficient is 
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As we have mentioned, the continuous models (Langevin equation and Fokker- 
Planck equation) can be regarded as infinite mechanochemical states cases of the 
discrete model, i.e. the periodic one dimensional hopping model. Therefore, the 
mean velocity and effective diffusion coefficient of continuous cases can be obtained 
directly by calculating the limits of which in discrete case ([6]). In the following, we will 
obtain the limits Vd and of Vn and -Dat by formal calculation. Theoretically, Vo 
and Vl, Vp are equivalent. Though Dp, and Dp, Dp are not theoretically equivalent, 
since Dn,Dp and Dp are all obtained under some assumptions, numerical results 
indicate that the differences among them are very small. 



2 The limit of velocity of the hopping model 



Firstly, we consider the large mechanochemical state N limit of the velocity of 
one dimensional hopping model. 
Due to the detailed balance 
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where Xi = 
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In the large N limit, Uj ~ (x)^ ~ Dj^/{Axy (for detailed discussion, see [21} [6]). 
Therefore 
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Thanks to the periodicity of the transition rates Uj and Wj ([2 

^ ^ ^ Wi+i [ (j){xN) - 4>{xo] 

= exp ( ^ — -j-^ . 
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(fig]) ([2U]), it is easy to obtain 
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where D = hnijv^oo -Dat. One can easily verify that, if D = D then Vd 
(see dHD (till))- In the following, we always assume that D = D. 
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(21) 



Vl = Vf 



3 The limit of effective diffusion coefficient of the 



hopping model 



Using the similar discussion as in the above section, it can be easily verified that 
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In the large biochemical state number N limit. 
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For the sake of the simphcity, we define 
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It can be readily verified that 



1 /"-^ 

Rn ~ =^ I f{x)dx 

N-l 



Moreover 



so 



UjTjSj 



j=0 

N-l 



S[Xj) fti = 



g{xj 



DAI) 

4-n J-iJ-'N Jo 



i=o 

iV-l 



+ l)rfe+,+i ^J2^k + l)f{xk+j+i) - 

k=o k=o 

= V Xk+ifixj + Xk+i) =- ^ -=- / zf{xj + z)dz 

t._n J-'J-'N Jo 



N-l N-l 



j=0 k=0 



N-l 
3=0 



Dn^ 



N 



LD 



N Jo 



zf{xj + z)dz 



^ N-l .L 

^^gixj) I zf{xj + z)dz 

D.T .-n Jo 



'N j=0 

N ''^ 



LDj^ Jo 



L 

g{x) I zf{x + z)dz 





dx 



(25) 



(26) 



(27) 



(28) 



7 



Substituting dZID dSSD dZHD ([SHD into (Q, one obtains 
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From (fT2l [T3|) and (!24l) , one can know that 
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Substituting fl3T|) into fl30l) . we finally get the limit of the effective diffusion coefficient 
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By the way, the hmit of the randomness parameter RAND^ = is 
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4 Numerical results and discussion 

As we have known that the hmit Vo of velocity of the hopping model is the same 
as Vl and Vp. However, it is easy to understand that D^, would be different from 
Dl,Df because they are obtained under different assumptions. 

To indicate the accuracy of the formulation (|32|) . we present the numerical results 
of the formulations ([H]) ( fT5l) and ( l32l) in Figure [1] and [21 From the numerical results 
we can find that the formulation ([32D is accurate enough, especially for the cases 
with small external force F^xt- Which implies that the methods used to derived the 
formulation fl32l) is reasonable. From the derivation we also can know more about the 
relationship between discrete models and continuous models of stochastic motion of 
microscopic particles. 

In conclusion, we have provided an analytical formulation of effective diffusion co- 
efficient of the stochastic motion of microscopic particles. The numerical comparison 
with the formulation obtained in the framework of the overdamped Langevin dy- 
namics and Fokker-Planck equation indicates that our analytical formulation is very 
accurate. Moreover, the methods used in this research can be further used to get 
more results about the stochastic motion. Through the discussion in this research, 
the relationship between continuous models and discrete models has also been made 
clear. 

Acknowledgments This work was funded by National Natural Science Foundation 
of China (Grant No. 10701029). 



9 



Figure 1: Left: figures of effective diffusion coefficient as a function of free diffusion 
coefficient D := ksT/C,, wliere tfie red line is obtained by tlie formulation f|T5l) . tlie 
green stars are obtained by the formulation ([8]) and the blue circles are obtained by 
the formulation fl52]) . Right: the difference between the formulations fl52]) ([HD and 
( IT5|) . where Dp denote the results of formulation (ITS!) . Dl denote the results of the 
formulation ([8]) and denote the results of the formulation (132!) . In the simulation, 
the potential ^{x) = Uosin{27!-x / L)) — F^xtX and ^ = Uq = L = 1 (see ^IH]), F^xt = 0.5. 
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Figure 2: Left: the relationship between the effective diffusion coefficient and the 
external force Fext- Right: the differences between the formulations (IHl) ( |T5l) and 
(!32|) . In the simulation, the potential = Uosin{2TTx/L)) — F^^tX and C, = Uq = 
L = 1 (see [18]). 
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